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We study cosmological perturbations around self-accelerating solutions to two extensions of non- 
linear massive gravity; the quasi-dilaton theory and the mass- varying theory. We examine stability 
(— ^ ■ of the cosmological solutions, and the extent to which the vanishing of the kinetic terms for scalar 

0~4 \ and vector perturbations of self-accelerating solutions in massive gravity is generic when the theory 

5— ( ' is extended. We find that these kinetic terms are in general non- vanishing in both extensions, though 

^^ there are constraints on the parameters and background evolution from demanding that they have 
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the correct sign. In particular, the self-accelerating solutions of the quasi-dilaton theory are always 
unstable to scalar perturbations with wavelength shorter than the Hubble length. 



5^, I. INTRODUCTION 

Recent years have seen the development of a non-linear theory propagating the five degrees of freedom of a massive 
" ' • graviton (dRGT theory [1, 2], see [3] for a review), without the Boulware-Deser ghost [4, 5]. This theory admits 
(^s^ . self-accelerating solutions [6-12], in which the universe is de Sitter without a cosmological constant in the action. 
'^ \ The Hubble scale of these self-accelerating solutions is of order the mass of the graviton. Having a light graviton is 
(<— ^ I technically natural [13, 14], so these solutions are of great interest to account for cosmic acceleration in the late-time 
_il ' universe. 

f^ ■ Given any non-trivial solution, it is natural to ask how the perturbations around it behave, and in particular 

whether there are interesting new effects in the propagation of the associated degrees of freedom. The perturbation 
theory for the self-accelerating solutions of dRGT has been studied in [15-19]. Freedom from the Boulware-Deser 
ghost means that around any background, at most five degrees of freedom propagate. Around a homogeneous and 
isotropic cosmology, these take the form of one transverse-traceless tensor, one transverse vector and one scalar. Even 
d • though the Boulware-Deser ghost is absent, the kinetic terms of these degrees of freedom can potentially have the 
wrong sign, in which case they are ghosts, signaling that that particular background is unstable. 

In fact, around the self- accelerating solutions of dRGT theory, the scalar and vector degrees of freedom have 
vanishing kinetic terms [15, 20]. This result could have critical implications for the cosmology, and it is important to 
understand the extent to which it is a generic result in these types of models. There are several avenues one might 
consider. Quantum mechanically, kinetic terms may be generated by loops. Determining the sign of such terms, and 
hence whether the background propagates ghosts, is then a difficult question whose answer depends in general on 
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the details of the matter propagating in the loops. If we wish to restore the kinetic terms at the classical level, one 
avenue is to move away from homogeneous and isotropic cosmologies [21, 22]. Another is to keep homogeneity and 
isotropy, but change the theory by adding more degrees of freedom. In this paper, we take the latter approach. We 
study cosmological perturbations in two extensions of dRGT theory: quasi-dilaton massive gravity, and mass- varying 
massive gravity. 

Massive gravity admits an extension to a theory with a global scale symmetry through the inclusion of a specific 
scalar field, dubbed the quasi dilaton [23]. We examine self-accelerating background solutions to the quasi-dilaton 
model, and consider the behavior of cosmological perturbations (see also [24]). Wc find the conditions under which 
the backgrounds are free of ghosts. We find that there is always a ghost-like instability in the scalar sector, for 
fluctuations of physical wavelength shorter than the Hubble radius. 

The other extension we consider is mass-varying massive gravity, the theory obtained by promoting the mass to a 
scalar field [10, 29]. We study the behavior of cosmological perturbations in this model, and find the conditions under 
which the backgrounds are free of ghosts. 

II. QUASI-DILATON THEORY 

We start with the quasi-dilaton theory, introduced in [23] . The action governs a dynamical metric g/^^i, and a scalar 
field a,^ 
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The part of the action which provides the mass to the graviton is 

/:3 = ^ {[K-f - 3 [K] [K.^] + 2 [/C3] 
A = ^([/C]*-6[/C]^[/C2]+3[/C2]^ 
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where square brackets denote a trace. While these expressions are similar in form to the dRGT theory, in the case of 
the quasi-dilaton theory, the building block tensor K is defined as 
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where g^u is a non-dynamical fiducial metric. The theory is invariant under a global dilation of the space-time 
coordinates, accompanied by a shift of a. This symmetry rules out a non-trivial potential for a. 
Throughout the analysis, we choose the fiducial metric to be Minkowski, 



g^vdx^ dx'^ 



dt^ + Sijdx"- dx^ 



(4) 



A. Background equations of motion 



For the physical background metric, we adopt the flat FRW ansatz 

g^ydx^" dx" = -N{t)^ dt^ + a{tf ^^Jdx^ dx^ 



(5) 



^ We thank the authors of [23] for pointing out the existence of the ^ term in the most general quasi-dilaton action. 



To obtain the background equations of motion, it is convenient to introduce time reparametrization invarianee, so 
that we may write a mini-superspace action. We replace the fiducial metric with 



g^,ydx^' dx" = - f'{tfdf + S^jdx' dx^ , 

where f{t) is the Stiickelberg scalar [13], and unitary gauge corresponds to the choice f{t) ~ t. 
The mini-superspace action is 
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+f'a^X [{X - 1) (3 - 3(X - l)a3 + {X - ifa^) +iX^] 
where V is the comoving volume and we have defined 
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oCr/Mp 



In addition, to simplify expressions later on, we define 

Varying with respect to / and then choosing unitary gauge, we obtain a constraint equation: 
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The Friedmann equation is obtained by varying with respect to the lapse TV, 
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and varying with respect to the scale factor, a, yields the acceleration equation. This may be combined in a linear 
combination with (11)) to yield the simpler equation 
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Finally, the equation of motion for a is 
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(13) 



Since time reparametrization invarianee was introduced with the Stiickelberg field /, there is a Bianchi identity 
which relates the four equations. 



SS . SS ■ ,dSS . SS ^ 



(14) 



Therefore one equation is redundant with the others and may be dropped. In discussing solutions we will generally 
choose to drop the acceleration equation (12), although we will use it to simplify expressions for the perturbations. 



B. Self-accelerating background solutions 

We now discuss solutions, starting with the Stiickelberg constraint (10). Integrating this equation gives 

x|(l -X) [3 - 3{X - l)a3 + {X - ifai] - ^ X^\ = — x constant. (15) 

In an expanding universe, the right hand side of the above equation decays as a^'*. Thus after a sufficiently long 
time, X saturates to a constant value Xsa, corresponding to a zero of the left hand side of (15). These constant X 
solutions lead to an effective energy density which acts like a cosmological constant. As pointed out in [23], there are 
four such solutions for which X is constant. Of these, X = implies a — >■ — oo, and as in [23], we drop this solution 
to avoid strong coupling.^ What remains are the three solutions to the cubic equation ■^ 

(1-X) [3- 3(X- 1)^3 + (X- 1)2^4] -^^3 =0. (17) 

For these solutions, we write the Friedmann equation (11) as 

(3-|)i/'=A + AsA, (18) 

where the effective cosmological constant from the mass term is 

Asa = ml (Asa - 1) [-3Asa + 6 + (Asa - 4)(Asa ~ 1)^3 + (^sa - 1)'«4] ■ (19) 

The Friedmann equation (18) also provides a condition on the parameter uj; on the self-accelerating solutions, one 
needs to have w < 6 in order to keep the left hand side of the Friedmann equation (18) positive. This ensures that 
when ordinary matter is added to the right hand side, we will have standard cosmology during matter domination. 
(Although we do not include matter fields, the sign of the matter energy density can be determined by replacing the 
bare A with p/M^.) 

Finally, on the self-accelerating solutions, with constant H specified by (18), the equation of motion for a fixes the 
ratio r = a/N . From Eq.(13), we obtain 



c^g^(AsA-l 

SA - 1)2 - 2(A; 

Here, to simplify the expression we have used the Stiickelberg equation (15) to eliminate a4 



'^^ ^ + ml As\ [a3(AsA - 1)^ - 2(Asa - 1) - C^IJ ' ^^^^ 



C. Perturbations 

To find the action for quadratic perturbations, we expand the physical metric in small fluctuations Sg^^i, around a 
solution gjjy , 

9^^^ = 9^^} + %^ , (21) 



^ As we discuss at the end of Appendix A3, the remaining solutions also lead to strong coupling in the vector and scalar sectors, when 

5 = and the parameters 03 and 04 are such that X ~ 0. 
^ For the choice § = 0, the system simplifies as one of the solutions to Eq.(17) becomes XgA = 1, while the remaining two are 



I 3 03 + 2(14 ±1/9 a? - 12^4 
X± ^ XsaI ^, = -^ . (16) 

In this special setting, the solution X = 1 is uninteresting: the effective cosmological constant from the mass term is zero and in 
the present scenario, the background becomes equivalent to a de Sitter universe driven by a (bare) cosmological constant 6A/(6 — w). 
However, in the presence of matter fields and no bare cosmological constant, this solution asymptotically approaches a Minkowski 
background and is unstable [23]. 



and keep terms to quadratic order in 5g^^. 

We break the perturbations into standard scalar, transverse vector and transverse-traceless tensor parts, 



55oo = -2Af"$, (22) 

5goi 



= Na (Bf + d,B) , (23) 



/»r + ^(^»^J + ^J^n + 2 6^J * + (^d.d, - i% did'^ E 



(24) 



Sgij = a^ 

where 

d%f;^ = hr^ = 0, d'Bj = Q, d'Ej^Q. (25) 

We then introduce the perturbation of the scalar via 

(T = a'^^^ + Mp Sa . (26) 

We perform the entire analysis in unitary gauge, so that there are no issues of gauge invariance to worry about, and 
no need to form gauge invariant combinations. We write the actions expanded in Fourier plane waves, i.e. V^ — > — A:^, 
(fix — > (fik. Raising and lowering of the spatial indices on perturbations is always carried out by (5*^ and dij. 

D. Tensor perturbations 

We begin by considering tensor perturbations around the background (5), 

a'/^r ' (27) 



where d'^hf^ = hf^ * = 0. The tensor quadratic action reads 



S = ^ fd^ka^Ndt(^\hjn'-(^ + M'aw]\hJn']^ (28) 



where the mass of the tensor modes is given by 

To obtain this, we have first used the background acceleration equation to eliminate any terms with a. Then we 
have used the self-accelerating branch of (15) (at late times when the right hand side is zero), the Friedman equation 
(11) evaluated on the self-accelerating solution (i.e. X = 0), and the a equation (13) evaluated on the self-accelerating 
solution (i.e. X ^ H = 0), to eliminate A, 0:3 and a4. 

The tensor mode always has correct sign kinetic and gradient terms. However, it will be tachyonic if the mass 
term is negative: M^.^ < 0. The stability of long wavelength gravitational waves is thus ensured by the condition 
MQyy > 0. Nevertheless, even if this condition is violated, the tachyonic mass is generically of order Hubble, so the 
instability would take the age of the universe to develop. 

E. Vector perturbations 

We next turn to vector perturbations. 
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Sg„, = Na Bt , %, = - {d,E^ + d,Et ) , (30) 



where d'^Bf = d^Ef = 0. The field Bf enters the action without time derivatives, so we may ehminate it as an 
auxiliary field using its own equation of motion (again we are using the equations of the background self-accelerating 
solution to eliminate A, a^ and a^) 

4wa2iJ2 + 2/c2(r2_l) N ' ^ ' 

Once this is inserted back into the action, what remains is a system of a single propagating vector, 

^_^ [ ^■iu^'i AT ^^ C^V \i^T\2 k^'hw \T^T\2 

where 



S=^ I d'ka''Ndt[^^\Et\^~'^^f^\Et\^] ^ (32) 



and Mgw is the mass of the tensor modes as in (29). 

From Eq.(33), we see that for (r^ — l)/w < 0, there exists a critical momentum scale, kc = aH J -^^^ , above which 
the vector becomes a ghost. In the case (r^ — l)/a; > 0, the kinetic terms of vector always has correct sign and thus 
there is no such critical momentum scale. In the first case, stability of the system requires that the physical critical 
momentum scale, kc/a^ be above the ultraviolet cutoff scale of the effective field theory, h-uv i i-e. 
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^Iv < -, ^ - if ir' - 1)/^ < . (34) 



To determine whether the vector modes suffer from other instabilities, we define canonically normalized fields, 

Sr^^TvEf, (35) 

in terms of which, the action (32) reads 

S=\jd'ka'Ndt (^-l|ff|2_^2^|£T|2^ ^ (36) 



where the dispersion relation of the modes is given by 

^^ = (1 + <f) M'aw - " ;, '::,,7 ' > (37) 

and we have defined the dimensionless quantity 



g 2g2(i^4g 2) 

(l + g2)2 



The second term in the dispersion relation (37), which originates from the time derivatives of 7Vj is always of order 
0{H'^), provided that q^ > 0. Therefore, in this regime, this term does not introduce instabilities faster than the 
Hubble expansion rate. Moreover, if q^ < 0, the no-ghost condition (34) imposes jg^j < (fc^/a^)/A^y. Thus, for any 
physical momenta sufficiently lower than the cutoff scale of the effective theory, the second term in (37) does not lead 
to any visible instability, i.e. the growth rate of any instability (if any exist) is at most of the cosmological scale. 

The vector modes may potentially suffer from a gradient instability arising from the first term in (37), if M^y^, < 
and g^ > 0. The growth rate of this instability can be made lower than or at most of the order of the cosmological 
scale for all physical momenta below the UV cut-off Auv , provided that 

A^v ^ ^^^ , if (r' - 1)/^ > and M^^ < . (39) 



F. Scalar perturbations 



Finally we consider the action quadratic in scalar perturbations, 
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The scalar sector should consist of two dynamical degrees of freedom: the scalar field and the longitudinal mode of 
the massive graviton. The perturbations $ and B stemming from Sga and Sgoo are free of time derivatives, and so 
we eliminate them as auxiliary fields using their equations of motion: 
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Inserting these back into the action, we obtain an action with three fields, ^P, E and 6<j. Since the "sixth" de- 
gree of freedom (which would come from the Boulware-Deser instability in generic massive theories) is removed by 
construction, there is another non-dynamical combination, which we determine to be 



^^—{^ + 5a) 



We also define an orthogonal combination. 



5a = 



V2fc^ 



{-^ -5a) 



(44) 
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With these field redefinitions, the action can be written in terms of ^, 5a and E, with no time derivatives on ^. The 
latter is therefore auxiliary and can be eliminated via its equation: 



/ , 24a2i72 2a2ij2fc2 r[48_(g_^)^]^_^2|\ 
* = -fc^ -, H ^-n^ 7^ -— ^^ 4r4 ^ 5a - 



2V2k^E 



~2a^H 



r(r-l) [4A:2_ (6-w)a;a2iJ2] (r_i) 

3 {6~uj)[2k^{r-l) + 3uja^H^]\ 5'a 



3 [4 P _ (6 - w)a; a2 iJ2] 
V2{6~uj)Pa^H E 



(46) 



^r [4fc2- (6-w)wa2i/2](r- 1) y iV 3 [4 k^ - {6 - oj)uj a^ H^] N ' 
Using this solution in the action, and introducing the notation Y = {5a, E), the scalar action can then be written as 
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where A^ is a real anti-symmetric 2x2 matrix, and /C and fl^ are real symmetric 2x2 matrices. (Note that there is 
no loss of generality in taking A4 anti-symmetric, since the symmetric part can be absorbed into H,^ by adding total 
derivatives). For now, we focus on the kinetic terms. The components of the matrix K. are 
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FIG. 1: The stability of the scalar sector implied by the determinant of the kinetic matrix (49). For modes with k/{aH) below 
the solid line, the determinant is positive, so there no ghost degrees of freedom (see Eq.(Al) for the field basis in which this is 
manifest). On the other hand, above the solid line, one degree of freedom has a positive kinetic term while the other is a ghost. 



For the case at hand, it is sufficient to study the determinant of the kinetic matrix /C to determine the sign of the 
eigenvalues.^ The determinant takes the comparatively simple form. 



det/C 



3M4fc6a;2a4^4 



H^- 



4fe2 



(r-l)2 



(49) 



The sign of the determinant is determined by the sign of the quantity within the square brackets. First note that the 
determinant is always negative if u; < 0. Along with the condition for a realistic cosmology obtained from (18), the 
range of allowed a; is thus 



< w < 6, 
in agreement with [23]. In order to have no ghosts in the scalar sector, we need (See Figure 1) 

k ■\/a;(6 — uj) 



aH 



< 



(50) 
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Gencrically, the right hand side of the inequality (51) is of order 1. This implies that for modes with physical 
wavelengths that are smaller than cosmological scales, one of the two degrees of freedom is a ghost. In other words. 



parametrically, there is an instability in the scalar sector at physical momenta above H 



As shown in Ap- 



pendix A, both the physical momenta and the energies of those ghost modes near the threshold are not parametrically 



higher than H - 



and thus are below the UV cutoff scale of the effective field theory. This signals the presence of 



ghost instabilities in the regime of validity of the effective field theory. 

We end this section by comparing our results to those found in [23, 24]. Noting that at the level of the kinetic 
matrix (48) the only scale other than H is the momentum, the limit iJ — > is equivalent to considering modes with 
wavelengths much shorter than the Hubble radius, i.e. k ^ aH, which is in contradiction with the no-ghost condition 
(51). In this limit, the kinetic matrix then becomes 



/C 



M^ UJ k^ 




aH 



(52) 



* The absence of ghosts requires that both the determinant and the trace are positive. On the other hand, only one of these being negative 
is enough to deduce the existence of a ghost degree, which happens to be the case for the current system. For a detailed diagonalization 
treatment, we refer the reader to Appendix A. 



which has one positive and one zero eigenvalue, as in [23] (See Appendix B for a more detailed comparison). In 
other words, the apparent stability of the self-accelerating solution is due to the loss of the dynamics of the ghost 
degree of freedom in the short wave-length limit, and so the decoupling limit is not sufficient to determine stability, in 
agreement with [23]. In [24], only the super- horizon limit A; — > is considered, so the instability which appears only 
for physical wavelengths < Hubble is not visible in this limit. 

G. Higher derivative terms and UV sensitivity 

The quasi-dilaton theory is governed by the global scaling symmetry described in [23]. The action (1) includes all 
possible terms compatible with the symmetry, with up to two derivatives, and we found there was no way to render 
the scalar perturbations of the self-accelerating solutions stable at all momenta. 

However, beyond two derivative order there are many more terms compatible with the symmetry. These higher 
derivative terms can be thought of as encoding UV effects from whatever physics completes the theory. Among the 
possible higher-derivative terms, we will focus here on two classes of distinguished interaction terms which will not add 
new degrees of freedom. There are the Goldstone-like terms of the form ~ (da)"", and the three possible non-trivial 
covariantized Galileon terms, of the form ^ (9cr)^(9^o')" + ■ ■ ■ [25-28]. The strong coupling scale of the quasi-dilaton 
on flat space is A3 ^ (Alpmiy^^ [23], so it is natural for the Galileon-like terms to appear suppressed by this scale. 
The Goldstone-like terms should carry the scale^ A2 ~ (MprngY/^ :$> A3. One can repeat the calculation of the 
perturbations including these terms, in the hopes that the fluctuations can be stabilized at short scales k ^ H. 

The Friedmann equation now becomes: 

3-'^ + 3g3h')H^ + lml[f{h')-2h\nh^)]=A + AsA, h^^, (53) 

"^3 



2 ■"■■'"■ W ■2--9L-'v-' -.-.,,.- ,j - , -.^, --^ 



where 173 is the dimensionless coupling for the cubic covariant Galileon and we have chosen the form Aj f{x)/2 as the 
Goldstone-like term, where x = — (V<t)'^/A|. On the other hand, the constraint equation and the value of AgA remain 
the same. (We have omitted the quartic and quintic Galileon terms for simplicity.) There are still self-accelerating 
solutions with H ^ mg so the existence of these solutions appears insensitive to the UV effects encoded by the higher 
derivative operators. For a sensible cosmology, H^ determined by the Friedmann equation should be an increasing 
function of the bare cosniological constant A (which represents the matter energy density in our setup) . Demanding 
this, we obtain the condition 

6 - w -h 12<73/i2 - /'(/i2) _ 2h^f"(h^-) > . (54) 

The determinant of the kinetic matrix JC for scalar fluctuations changes by order one, 

detK, — 7 — — ^ , (55) 

but the determinant is still always negative for sufficiently large momenta, provided that the condition (54) is satisfied. 

The quartic and quintic covariant Galileon terms can render the determinant of the kinetic matrix for scalar 

perturbation positive for large momenta, depending on the values of the coupling constants. However, in this regime 

of parameters, the tensor and vector modes acquire negative kinetic terms. Moreover, after explicit diagonalization 



^ The reason the Goldstone-Uke terms carry a higher scale is because the A3 decoupling limit of the theory has an enhanced Galilean 
symmetry [23], which the Goldstone-type interactions are not invariant under. This means that they will not be generated in the 
decoupling limit, so whatever the quantum corrections to these operators are in the full theory, they should not survive in the decoupling 
limit, i.e. they should be suppressed by a scale higher than A3. 
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of the kinetic matrix one can show that there are two ghost modes in the scalar sector, provided that H^ , determined 
by the Friedmann equation, is an increasing function of A. 

Thus, the form of the dispersion relations for the perturbations depends on and receives order one correction due 
to UV effects, but the presence of the ghost seems to be a robust feature. 

III. VARYING MASS THEORY 



We now turn to the varying mass theory, obtained by introducing a scalar into dRGT theory and allowing the 
graviton mass to be a function of this scalar. This theory was first considered in [10], and further studied in [29]. 
The action is 
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a X v~5 



Mi 
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We have further generalized to allow the dRGT parameters a^ and 04 to depend on the scalar a. The part of the 
action which provides mass to the graviton takes the same form as in Eq.(2), but here the building block tensor JC is 
the same as in dRGT theory [2], i.e. 



ICt 



s^.-{vF^y 



(57) 



One of our goals is to compare our results with the analogous analysis of perturbations in the dRGT theory. Since 
the original theory does not allow flat solutions for a Minkowski reference metric, it is necessary to adopt a more 
general form. We therefore extend the fiducial metric to be an arbitrary spatially flat homogeneous and isotropic 
metric, 



g^j^^dx'^dx'' = -n{tYdr + a{tyS^jdx^dx^ 



(58) 



A. Cosmological Background Equations 

We first study the cosmological background equations (see [29-34] for more on background cosmological solutions 
to mass-varying massive gravity.) For the physical background metric, we adopt the flat FRW ansatz 



g^vdx^ dx^ 



N(tydr + a{ty6ijdx'dx^. 



(59) 



To write the mini-superspace action, we introduce time reparametrization via a Stiickelberg field /(t), by replacing 
the fiducial metric with 



g^^dxf'dx" ^ -n{f{t)Yf'{tYdt' + a{f{t)Y5ijdx'dx^ 



Unitary gauge corresponds to the choice f{t) = t. 
The mini-superspace action is 



V 



^ dt I 'SM, 



N 



+ m'[NF~fn{f)G 



where V is the comoving volume and 



-N-^a"^ - NV{(j) -N MlK 



"3, 



"4, 



F^a{a- a(/))(2a - a(/)) +^f{a- «(/))2(4a ~ «(/)) + -j{a - a{f))\ 
G ^ a\a - «(/)) + a^a{a - a{f)Y + ^(a - a{f))\ 



(60) 



(61) 

(62) 
(63) 
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In the following, for clarity, we will use the definitions 



JNa a na 



NX 



(64) 



and we will omit the dependence of the functions nig, a^, 04 and V on the field value a. (We also caution the reader 
that the above definitions of X and r are different than the ones in the quasi-dilaton theory, which we introduced in 
Section II.) 

The equation of motion for the temporal Stiickelberg field / is 



1 5S 



ZNQNn 6f 



f=t 



^i{->'('-'> 



1- (X - I) as + ^{X - 1)'^ ai 



a^HmlX 



3 - X(2 + r) + (X - 1)((1 + 2r)X - 3)^3 - {X - lf{rX - 1)0^4 



= 0. 

(65) 



The Friedmann equation is obtained by varying the action with respect to N , 



1 5S 



M^a^ SN 



3H^ -A 



f=t 



^ (2^ + ^) - "^' (^ - 1) [-3{X-2) + {X-4){X-l)as + {X~lfa,] ^0, 



(66) 

and by taking a variation with respect to a, we obtain the dynamical equation which, after forming a linear combination 
with (66), can be expressed as 



1 SS 



2,M?:Na? 5 a 



1 5S 



j^. My SN 



2H a^ 



f^, N M2 7V2 9 



{r~l)X\;3-2X + {X - 3){X - 1)0^3 + {X ~ ifa^] = . 

(67) 



Finally, the equation of motion for a is 



1 5S 



a^N 5 a 



f=t 



N dt \N N 



Ml ml {X - if I 4 (4 - X(l + 3 r)) + a'^{X - I) {r X ~ I) 



2 m' 



3(X(r + l)-2) 
X -1 



{X{1 + 3r) ~ 4)^3 + {X~ l){rX - l)a4 



0, (68) 



where a prime denotes differentiation with respect to a. 

It is convenient to cast these equations into a more familiar perfect fluid-like form, by defining the following 
quantities 

p„ = M^m] {X - 1) [-i{X -2) + {X-A)(X- 1)^3 + {X - ifa^] , 

p„ = M^ml [6 - 3X(r + 2) + X^{1 + 2r) - {X - 1)(4 - X(3r + 2) + rX^)a:i - {X - lf{rX - 1)04] , 



a 

N 



Q = Mlml^{X-lf{a',{4~X{l + 3r)) + a'^{X-l){rX-l) 



2 m' 



3{X{r+l)-2) 
X -1 



- (X(l + 3r) - 4)a3 + (X - l)(rX - 1)04 



Pa 



Pa 



2iV2 



+ V, 



2iV2 



V, 



(69) 
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in terms of which Eqs. (65)-(68) can be re-written, respectively, as 



Pm 

N 



2,H{p,n +Pm) = -Q, 

P 



H 

N 

Pa_ 

N 



1 



2M2 



[{Pa +Pa) + [Prn + Pni)] 



3H{p,+p,) = Q. 



(70) 



From here on, we will use these forms of the cosmological background equations. 

B. Tensor perturbations 

We proceed with the perturbation theory in the same manner as described in Section II C. We start with tensor 
perturbations around the background (59), 

(71) 







5g.,-a^hJ^, 


where d^hf/" = hf^ ' ^ 0. 






The action for the tensor perturbations reads 






Af2 r 


( 1 |/,TT|2 


-(5 + a4vv)|/^ST)> 


where the mass term is 






^^GW ^x - 1)2 


mliX 


1) ^™ 1 


/ 1 2X \ p™+p™ 


\r-l X-l) Af2 



(72) 



(73) 



To obtain this, we have used the background acceleration equation (67). In the case of self accelerating solutions 
[11, 15] of the standard dRGT theory, i.e. when (7 = and pm — —Pm = const., the last term in (73) drops out of the 
calculation, and the mass reduces to the one found in [15]. We also stress that Mqw here is different than the one 
defined for the quasi-dilaton theory in Section II. 

C. Vector perturbations 

Next we consider transverse vector perturbations to the metric 



5go; = NaB. 



T 

i ' 



6g., = -{d.Ef + d,Ef)., 



(74) 



where d^Bj = d'^Ef ~ 0. The field Bi appears without time derivatives and may be eliminated as an auxiliary field 
using its own equation of motion. 



B 







a 


2 


1- 


' fc^A/i^°r-2-l)(Pm+P™)J 



N 



(75) 



Once this solution is inserted back into the action, what remains is an action for one dynamical vector 

S^^A jd^ka'Ndt {l^^\Ej\--^^^\EjA , (76) 
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where 

1,2 / y2 (^2 _^-^ j^f2\ -1 

(77) 

and the tensor mode mass Mqw is as defined as in (73). For the self-aecelerating solutions of the dRGT theory, where 
(7 = and p„i + Pm = 0, the vector kinetic term vanishes, in agreement with the results of [15]. 

From Eq.(77), we see that for (r^ — l)/{p„i +Pm) > there is a critical momentum scale above which the vector 
modes become ghosts. On the other hand, in the opposite case with (r^ — l)/{pm +Pm) < 0, the timelike kinetic term 
of vector modes always has the correct sign, and thus there is no such critical momentum scale. The stability of the 
system thus requires that such a critical momentum scale should be either absent or above Ajjv j the UV cutoff scale 
of the effective field theory, i.e. 

In order to determine the stability conditions and the time scale of potential instabilities, it is useful to use canonical 
normalization. However, the existence of several unknown functions and the lack of a simple background prevent us 
from performing the stability analysis in a complete way. On the other hand, assuming that the tensor modes have 
positive squared- mass (73), i.e. Mq^^ > and the vector sector is free of ghosts (78), we can still obtain sufficient 
(but not necessary) conditions to ensure the stability of the modes. These conditions cause the vector modes to damp 
with time, so tachyon-like instabilities can be avoided. It is convenient to perform a time reparametrization choosing 
the lapse function to be iV = a'^TV 

5U=a3r. = ^ Jd'kdt i^\Er\' - ^a^TvM'a^lEn^^ . (79) 

Hence, provided that the conditions MQ^, > and TV > are already imposed, the amplitudes of the variables Ef 
decrease as the universe expands if 

dt [a^TvM^^] > . (80) 

Demanding that this condition holds for all momenta below Aj/y, we obtain 

Afjyil-r') 3 

where 

A^l + ^—-\n(^^), B=l + ^—-\n(Ae^) . (82) 

8NH dt V ^ - 1 / 6iViJ dt ^ ^^' ^ ' 

D. Scalar perturbations 

We now move on to the scalar perturbations. In the absence of matter, we expect the sector to contain two degrees 
of freedom. The scalar parts of the metric perturbations are 



6goa = -2N^<^, 5go^ = Na^^B, 6g,j = a'^ 



2 5,,*+ { d.,dj - ^S^j did' ] E 



(83) 



while the scalar field is expanded as, 

a = cr'-^^ + Alp 5(7 . (84) 
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The perturbations $ and B coming from 5gai and (5goo carry no time derivatives and are non-dynamical, so we may 
determine them using their own equations of motion: 



B = -- 



Ml {r + 1) 



4 M4 m ^ (r2 - 1) + (/,,„ + p„,)ip, +p,-6 A/2 i/2) 
^'^"■"^^ 2MlH^E+{p,+p,-QMlH^)^ 



N 



-2H 



2 M2 ^ (r - 1) - 3 (p„, + p,n) * + 2 (r - l)(p, + p,) - 



* 2MpHir~l){p,+p,) 



N 



S& 



MpHN f r-1 \ {r-l){p,+p,-6MlH^)& 

' 6(H - HrX){p,n+Prn) -— (pa - Pa + 2 p,,^ ^ ' 



N 



MpN 



Sa}, 



$ = - 



M; 



#2 



4 Mp4 i/2 ^ (^2 _ 1) + (p^^ + p^^J(p^ +Pa-6 Af 2 i/2) 



3a2 
-2H 



ip,n+Prn)E-2MlH{r^-l) 



2M^^{r'^l)-3{p,n+Prn) 



E 

N 



-'?n I Frn 



Af 2 (r 



2M^^ir-l)-3{p,n+Prn) 



* 



* {Pni + Pm){Pa + Pa) 



{Pra +Pm)N 



2A/2(r-l)tT 



TV 
6(iI-iJrX)(p„,+p„0 



Mpcr 



S& 



-—ip^-p^ + 2p„,) -^ 



6a 



(85) 



Inserting these back into the action, we end up with a system of three degrees of freedom, ^P, E and Sa. Since the 
would-be Boulware-Deser ghost is removed by construction, there is another non-dynamical combination, which is 
found to be 



~ 1 / MpHN^ 
* = ^ * H ^ Sa 

\/2 V ^ 



We also define an orthogonal combination. 



Sa 



\/2fc2 



* ^ Sa 



(86) 



(87) 



The action can now be written in terms of ^, Sa and E, with no time derivatives on ^. The latter is auxiliary and 
can be eliminated with its own equation of motion. Thus, we obtain an action in terms of Sa and E of the form 



5': 



d^k 



■CL^ Ndt 



Y^ Y Y^ Y 

N N N N 



where Y = {Sa, E'), A^ is a real 2x2 matrix, and K, and fi^ are real symmetric 2x2 matrices. Note that by adding 
boundary terms, the mixing matrix M between fields and derivatives can be made anti-symmetric. 

The full kinetic matrix is rather lengthy, and so we will not display the full expression, other than to note that for 
the dRGT theory, where a ~ and pm = —Pm, it can be checked that /C = 0, consistent with the results of [15]. 
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Specializing to a Minkowski reference metric (H = H = 0) brings the kinetic matrix to a more manageable form: 



K-ii = 



(r-l)2 



18 {pm +Pm) 



2M^B + ?,a-^{p^+p,r.) 



iM^P 



2 M2 f^ (r - 1) - 3 {p,n + Pm) (Pa + Pa) 



JCu 



\/2A/2fc4 



[2 M2 fc2 + 3 a^p^ + p„)] ^4 M4 H^ ^ - (p,„ + p„)(p, + p, - 6 M2 i72) 

9r{pr,i +Prn) 



3(r-l) [2M2^ + 3(p™+p„0 



2Af2fc2 



2 M2 ^ (r - 1) - 3 (p™ + Pm) (Pa + Pa) 



IC22 = 



18 



[2 M2 fc2 + 3 a2(^^^^ + p^^)] ^4 M4 i/2 ^ _ (^„^ + p^„)(p^ + p, - 6 M2 ^2) 

9(pm +Pm) 



2Af2^ + 3(p™+p™) 



4M4fc4(p,+p,) 



with determinant 



det[/C] = 



a2[2 Mp2 ^2 + 3 „2(p^^ + p^^^)] ^4 M4 iJ2 ^ _ (p„^ + p,^J(p^ + p. - 6 M2 i/2) 
3 Af 2 q2 fc6 (^^ ^ p^)2 (^^ + p^ _ 6 M2 g2) 



(r-l)2 



4 M4 iJ2 f^ - (p„ + p„)(p, + p, - 6 A/2 H^] 



(89) 



(90) 



We stress that we have not specified any background solution up to this point; the only choice we have made is to fix 
the fiducial metric to be Minkowski. By requiring that the determinant is positive, we infer that in order to avoid a 
ghost degree of freedom, the momentum should satisfy® 



Pa+Pa 3\ lf_ Pin+Pm 



4 A/2 H'2 2 



Af2 



(91) 



Note that this condition should be imposed for all k in the regime < k/a < Ajjv, where Ajjv is the UV cutoff scale 
of the theory. 

In a regime in which we have a de Sitter like expansion, i.e. \H\ <C H^, this condition becomes even simpler, 

4 fc2 



R 



R-6 H'^a? 



>0 



(92) 



where R is defined in (78). Demanding that the condition (92) holds for all physical momenta k/a < Ajjv and 
supposing that Auy/H > 3/2, we obtain the no-ghost condition for scalar perturbations in the regime \H\ <C H^ as 



R>6 



(93) 



E. Consistency of stability conditions 

We now discuss the regions of parameter space in which the stability requirements we obtained in Eqs.(73), (78), 
(81) and (93) can be satisfied. The summary of the conditions is: 



Note that absence of ghosts requires that both the determinant and the trace are positive. However, for the scenario at hand, applying 
the field redefinitions given in Appendix C shows that one of the degrees of freedom always has a positive kinetic term. Thus, Eq.(90) 
is enough to ensure a healthy kinetic action. 
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i.) To avoid a tachyonic instability in the tensor sector, we need (from Eq.(73)) 

M^w>0- (94) 

a.) To avoid a ghost instabihty in the vector sector, from Eq.(78), 

in.) To avoid the unchecked growth of vector perturbations, from Eq.(81), 

iv.) To avoid a ghost instabihty in the scalar sector, from Eq.(93), 

i? > 6 . (97) 

Here, we have defined 

(p™+P™) ^^,^^1^^(RMgw)^ B=l + ^^ln(M'aw), (98) 

u2 A r2 ' ■' o AT u J+ \ „2 i / ' a at zj ,u \ uw J i \ / 



i/2 M2 ' 8NH dt V ^^ - 1 / ' 6A^iJ dt 

and have assumed that the UV cutoff scale Kuv is higher than 3/2 in units of H and that the expansion is de 
Sitter-like, i.e. \H\ ^ H^ (relevant for the scalar sector no-ghost condition). 

Note that if we satisfy the condition iv., then the condition ii. is trivially satisfied if r^ > 1, and that the condition 
Hi. is also trivially satisfied in this case if both A and B are positive. In more general cases, the above set of stability 
conditions is less trivial, but in principle there are regimes in which all of them are simultaneously satisfied. 

IV. DISCUSSION 

If the cosmologies of any of the recently proposed variations of massive gravity are to be of phenomenological use, 
it is crucial to understand the extent to which the theories propagate well-behaved, ghost free perturbations around 
their cosmological backgrounds. In this paper we have carried out this calculation for the cases of the quasi-dilaton 
theory and for the mass varying massive gravity theory. We find a host of constraints on these theories, primarily 
stemming from the requirement that ghost degrees of freedom not appear in the regime of applicability of the effective 
field theory. In the case of the quasi-dilaton theory, it can be seen that the stability found in the decoupling limit is 
an artifact of that particular limit and that, in fact, a ghost degree of freedom remains in the full theory. 
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Appendix A: Dispersion relation for scalar perturbations in the quasi-dilaton theory 



In this Appendix, we provide the details of the diagonalization procedure for the quadratic action for scalar per- 
turbations in the quasi-dilaton theory, starting from (47). 

1. Canonical normalization 



We first introduce a new field basis 



^1 



3^/uJ 



a^H'^ 



a; (6 — w) 



^1/2 



E + qV2 I 1 



2 (r - 1) fc2 



5a 



k^ M„ 



V6 



^ (E + 6V2Sa) , 



(Al) 



in terms of which the kinetic matrix (48) becomes diagonal and canonically normalized, so that the action is formally 

(fk 



S 



■a^Ndt 



N N N N 



with a canonical form for the kinetic matrix: 



/C 



Sign(l-A:2) o 
1 



where we have introduced the dimensionless (and time dependent) momentum via 

~ 2k 



a H y^uj (6 — uj) 
By adding appropriate total derivatives, the mixing matrix A4 can be made antisymmetric, reading 



M 



U kJl-k^ 

l--H(2r-l) ^ . 

6 ^ ^ l-fc2 



1 
-1 



Finally, the components of the symmetric matrix fl^ are 



(f^')ii 



(^(2)12 



l-P 



\r-uj-2){Sr-uj-S) + 



4{8r-uj-6) 
l-fc2 



12 



(l-fc2)2 



l--ksJ\l-k^\{M^awir-^) 



-H' 



2r(8r-w-10)+cj + 4- 



(2r-l)(8r-w-7) 3(2r-l) 



{n'h 



M'aw+^^^\M'awir-l) + 2H' 



(l-fc2) 

(2r-l)2^ 



(l-fc2) 



2\2 



<^^-^)-S+'- l_fc2 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



where MQ^, was defined in Eq.(29). For fc > 1, the mode Zi becomes a ghost, while for momenta fc < 1, both degrees 
of freedom are well-behaved. 

At low momenta (fc ^ 1), the action effectively becomes diagonal. In this long wavelength regime, both degrees of 
freedom have positive kinetic terms, while Mij = and the two eigenfrequencies are uif = and a;| = M^yy [24]. On 
the other hand, the ghost degree of freedom appears at momenta fc ^ 1, and so we must still diagonalize the system 
to determine the amplitudes of the frequencies in this regime. 
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2. Diagonalization 



To find the eigenfrequencies of tlie system, first note tliat since tlie matrices M and 51^ are time dependent, it is 
not possible to diagonalize tlie system at tlie level of the Lagrangian. On the other hand, the Hamiltonian can be 
written as a sum of decoupled oscillators, as shown in Ref. [35] and in the presence of ghosts, in Appendix D of Ref. 
[22]. We consider the two cases k <1 and fc > 1 separately. 

a. No ghost: fc < 1 



For the first case (denoted by subscript <), the kinetic matrix is unity and there is no ghost. In this case, we can 
introduce a rotated basis W = i?< Z, where the S0(2) rotation i?< satisfies 



This rotation allows us to remove the mixing, and the Lagrangian becomes 



(A7) 



N N 



.<.v.<^..<.. . (A8) 

As shown in Appendix D of Ref. [22] , the eigenvalues of the matrix 

nl = nl+MlM<, (A9) 

correspond to the actual eigenfrequencies. The matrix J7^ can be diagonalized by performing an S0(2) rotation ^<, 



f < ^< ^< = ^< (diagonal) 



where 



e< = 



and 

sin(2 6l<) 



(fi2 



<n2 



cos(0<) sin(0<) 
— sin(0<) cos(6'<) 



cos(2 6'<) 



(^1)22 -(^1) 



<;il 



The eigenvalues are then 

1 



{n^ 



<n2 



(AlO) 



(All) 



(A12) 



in%)n - {ni)22 



{n^ 



<n2 






{n^ 



<;ii 



(f^l)22 



1 



<J2 



2 \ , ^o2 \ , . / /o2 ^ t^2 



{ni) 



<;i2 



(r!^<)ii + (r!^<)22 + \ {nl),i - {Qlh2 



(^'<)l2 



(A13) 



b. One ghost: fe > 1 

In this regime (denoted by subscript >), the kinetic matrix has Lorentzian signature and the first mode is a ghost. 
Again, we introduce a rotated basis W = i?> Z, where the S0(1,1) rotation i?> satisfies. 



R> =NR>M>T], 



(A14) 
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with 77 = diag(— 1, 1). This rotation removes the mixing, and the Lagrangian becomes 



£>- ^ 77^ .. ..> V"> ' — > 



W'^ Rl {nl.+Ml.riM>)R>W 
In this case, the matrix we need to diagonahze is [22], 



which gives the actual eigenfrequencies 



For an SO (1,1) rotation given by 



^l = nl+Mlr^M> 



£,y Vti, ^> = 77 cj?. (diagonal) . 



e> 



where 

sinh(2 6'>) = 

the eigenvalues are 



2(f^^)l2 



cosh(0>) sinli(0>) 
sinh(6'>) cosh(0>) 



cosh(2 6l>) 



{nl),, + {nl) 



>;22 



[hi) 



2 

>;i2 



(A15) 



(A16) 



(A17) 



(A18) 



{^l),^ + {nl)2 



{nl)^i + {^i)22 



{n^ 



>)12 



(A19) 



(^>)i 



(^>)2 



(l^i.)22 - {ni)ii - \ {^lU + (^2^)22 



(^'>)l2 



{ni)22-{ni)ii + \ {nl)^i + {nl) 



j22 



-4 



(^'>)l2 



(A20) 



c. Comhinmg the two regimes 

Now that we have the necessary tools to diagonalize the system for the two regimes of momenta, we unify the two 
results. We first note that from Eqs.(A5) and (A6), we have 



which imply 



(n 



<)22 - -(f^>)22 , 



[inl),2]' = -mu- 



[(A<<)i2]' = -[(Al>)i2]', (A21) 



{nl)n = -(f^'>) 



>ni ' 



(f^<)22 



-(f^>)22, 



(f^<)l2 



(f^>)l2 



or 



(W<)? = (W>)?, (W<)2 = (^>)2 



(A22) 



(A23) 



Thus, it is straightforward to write down a unified expression for the dispersion relation, independent of the 
momentum regime of the modes. We obtain, 



P 



.12^-^<^H^ 



(8r2 + |-ll(2r-l)) +4 



(6-w)7-2 +Lu{r-l) + 3 

l-fc2 



9 



(l-fc2)2 



+ M2^(6-a.)(r-l) 



M. 



r^iV' 



2(6-w)fc2(i-_p) 



S2. 



(A24) 



20 



where a — (+) sign corresponds to the first (second) cigenmode, and we have defined 

6 + uj + 8{r - l)[r{6 - uj) - 3] 9 




B = H' 



j(6 - w) - 2 (r - 1)[4 (12 -uj)r-5uj- 12] + 



2(r-l)[6-w + 8(r-l)]-a; 



l-P (l-fc2)2 

r-w-7)(2r-l) 3(2r-l) 



l-fc2 



(l-fc2)2 



+ Af^^(r-1). 



(A25) 



3. Stability 



Since the modes with momenta fc > 1 are ghosts, we need to determine how serious this problem is. The ghost 
mode appears at (physical) momenta parametrically of the order of the Hubble rate. If the frequencies of these modes 
are larger than the UV cutoff of the theory, they are not within the regime of validity of the low energy effective 
theory and may be ignored. 

In the transition region where /c — > 1, the frequencies are 



3^2 



4(fc-l)2 

6 — w 



O 



1 



Mqw 



H 



M^yyir - 1) + —- [w(16r-5)+4r2 [2(6 - cj) r^ + 4 w r - 6a; - 9]] } + 0{k - I) {.A26) 
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We note that the problematic mode, right after k 
set of parameters 



1, has a very large frequency. As an example, we consider the 



A = 0, 



1. 



aa 



-10. 



Ui = 6 , 



^ = 0, 



mg<0.. 



branch , 



(A27) 



where "+ branch" corresponds to the positive sign solution in Eq.(16). These parameters lead to r ~ 1.01 and 
Mq^ ~ 0.71 |TOg|, which satisfy the stability conditions for the tensor and vector modes. For this example, we show 
the momentum dependence of the scalar dispersion relations in Fig. 2. As discussed in the paragraph after Eq.(A6), 
at low momenta, aj2 — >■ 0, while a;| — )■ M^-^. After the transition region, where aj2 exhibits divergent behavior, both 
modes increase with uP' o^ l? . The "light" mode, which becomes a ghost in the fc > 1 region, has (for this specific 
example) frequency uJi/ H oc 0(1), so that apart from in the immediate neighborhood of fc ~ 1, it cannot be integrated 
out from the low energy effective theory. 

Next, we consider a more general example, and extend our analysis to the a^, a^ parameter space. Instead of 
analyzing the immediate neighborhood of the critical point, we chose fc = 2, such that the frequency of the ghost 
mode becomes finite, while differing from the critical point value by an order one factor. In Fig. 3, we show the order 
of magnitude of tu'^ /H'^ for an example with A = 0, w = 1 and ^ = 0, in both branches defined in (16). In the regimes 
where wf > 0, the frequency is always of 0{H), and the ghost mode cannot be removed. On the other hand, we see 
that in a special region, the ratio \uJi\/H may exceed 102. The source of large jwij is related to a specific relation 
between parameters. 



a4 



-3(1 + 03) 



as > — 2 , + branch 
aa < — 2 , — branch 



(A28) 



which leads to X = and r oc X~2 _>. 00. Although we have excluded the line (A28) from our analysis, parameters 
close to this line lead to the large (negative) values we observe in the squared-frequency of the ghost mode. In 
principle, the value of \u!i\/H can be made arbitrarily large by tuning 0:3 and 04 to be close enough to this line. 
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FIG. 2: Plot of the dispersion relation of scalar modes versus the rescaled momenta k, for the example (A27), with A = 0, 
uj — 1, as = —10, a4 = 6, ^ = 0, rri^ < 0, in the positive branch defined in Eq.(16). The dashed line corresponds to Mq-^y /H^ 
which is the mass term for mode 2 (blue) in the low momentum regime. At the critical point k — 1, lj^ (red) diverges to — oo, 
and then becomes positive and finite after the transition. 



Mq^, 



Finally, we note that in general, close to line (A28), the parameter Mq^/ defined in Eq.(29) becomes negative, with 
-uj H^ + 0{X). In other words, even if the parameters are fine-tuned to remove the ghost mode in the scalar 
sector, the tensor and vector modes have a tachyonic instability (although its rate is at most of the Hubble scale). 
Additionally, the fact that X — > and r — > oo is also an indication of strong coupling, since the determinant of the 
kinetic matrix (49) is of the same small order of magnitude as the fine-tuning between 04 and 03 . 



Appendix B: Comparison of the scalar sector of the quasi-dilaton theory with results in the decoupling limit 



In this Appendix, we compare our results for the scalar sector of the quasi-dilaton theory with those obtained in 
[23] in a decoupling limit. We first identify our degrees of freedom with the ones used by Ref.[23]. By comparing 
perturbations of 5gij , wc find that 



* 



d'^b 



H 



3 nig a^ Trig 



An 



E 



(Bl) 



where ^ and E are the perturbations introduced by the decomposition (40) and Aq, b are the perturbations used in 
[23]. For the quasi-dilaton perturbations, 6a in (41) coincides with Q in [23]. 
We now turn to the action (47), which has the following kinetic term 



2 J N N' 

where the components of the kinetic matrix are given in Eq. (48) and the field basis is given by 

E J 

or, using Eq.(Bl), 



Y 



(B2) 



Y = 



(B3) 



/ 1 ( d^b 


-)] 


V2fc^ \3niga^ 


2b 


/ 



(B4) 
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, + Branch 






FIG. 3: Order of magnitude plots of the squared-frequency in Hubble units for the case A = 0, oj = 1 and ^ = 0, and for a 
benchmark momentum value k = 2. The left (right) column shows the positive (negative) branch of solutions given in Eq.(16), 
while the upper (lower) panel shows uji [002)- The color representation is as follows: Violet - (0,1] ; Blue - (1,10] ; Green - 
(10, 10^] ; Yellow - (10-, 10'^] ; Orange- (10'\ 10^] ; Red - (10*, 00) . The gridded region corresponds to oj^ < 0, while the thick 
black line corresponds to Q4 — —3(1 + Q3) [for the positive (negative) branch, 03 > —2 (03 < —2) only], where X ~ 0. 

with p = C + {H/m)Ao. 

Since the decoupling limit action in Rcf. [23] is given in the Z = {p, b) basis (up to the non-dynamical degree of 
freedom Aq), we transform our action via 



Y = RZ 




(B5) 



so that the kinetic term in the Z basis becomes 



ICz = R' ICR. 



(B6) 
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Before taking the decoupling limit, we consider the determinant of the kinetic matrix, given by 

2 



detICz = -- 



V2 



det/C . 



(B7) 



In other words, the momentum dependent ghost-free condition (51) is still valid. On the other hand, if we go to the 
decoupling limit, given by 



nig — > , 



the kinetic matrix in the Z basis becomes, 



/C, 



TT 

iJ — > , = finite , 

Too 



0(e) 0{e^) 



(B8) 



(B9) 



where e denotes the order of m^ and H. Thus, we see that at momenta comparable to and smaller than the expansion 
rate, one degree of freedom becomes a ghost, as we found in the main text. Therefore, this decoupling limit is not 
sufficient for determining the stability of one of the degrees of freedom. (Also see the discussion at the end of See. II F). 
This result coincides with the conclusion in [23], from the determinant in Eq.(B7). 

Appendix C: Diagonal Basis for the scalar sector of the varying mass theory 

In this Appendix, we diagonalize the kinetic matrix of the scalar sector in the varying mass gravity theory, studied 
in Section HID. Specifically, we want to show that the condition 

det/C >0, (CI) 

for the kinetic matrix is enough to ensure the absence of ghost degrees of freedom. 
The kinetic part of the action is given by 

J 2 N N 

where the components of the kinetic matrix are given in Eq.(89) and Y = (da, E). We now define a new basis, 

fc2 AL 



(C2) 






E + 6V2(l 3aMp™+P™; 



2fc2M2(r- 1) 



da 



V6 



^ E + 6 



V2Sa] 



after which, the kinetic terms become diagonal 



SB 



zl Zi 



^2 ^2 



with 



Kl 



— a Ndt \ Ki^ h K2 

2 \ N N N N 



Pm \ Pm 



"■ ^ \iMiH'i ~ 2) P 



K2 = l. 



(C3) 



(C4) 



(C5) 



Thus, the condition (91), obtained from the positivity of det/C, actually corresponds to the sign of the kinetic term 
of Zl, while Z2 always has positive kinetic term. 
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